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ABSTRACT 


During the course of this study, three tasks related to soil 
moisture sensing at microwave wavelengths have been undertaken: 

(1) Analysis of data at L, X and K 21 band wavelengths over bare 
and vegetated fields from the 1975 NASA sponsored flight experiment 
over Phoenix, Arizona, (2) modeling of vegetation canopy at micro** 
wave wavelengths taking into consideration both absorption and 
volume scattering effects, and (3) investigation of overall atmospheric 
effects at microwave wavelengths that can affect soil moisture 
retrieval. 

Data for both bare and vegetated fields are found to agree well 
with theoretical estimates. It is observed that the retrieval of 
surface and near surface soil moisture information is feasible through 
multi-spectral and multi -temporal analysis. It is also established 
that at long wa/elengths, which are optimal for surface sensing, 
atmospheric effects are generally minimal. At shorter wavelengths, 
which are optimal for atmospheric retrieval, the background surface 
properties are also established. 



1 . INTRODUCTION 


The feasibility of monitoring soil moisture utilizing passive 

microwave sensors has been tested through various truck -mount, airborne 

and space programs. At microwave frequencies the dielectric coefficient 

of dry soils is much less than that of water. As moisture is added 

to the soil, the dielectric coefficient increases. The thermal emissions 

of material are inversely related to their dielectric coefficients. Thus, 

the brightness temperature, defined as the apparent temperature measured 

by an antenna assuming that the emitting surface is that of a black body, 

is much lower for a moist soil than that of a dry soil at the same physical 

temperature. From a space platform the use of microwaves has the advantage 

over more conventional techniques that observations can be made at night 

and through cloud cover. Also, because of the relatively large skin 

depths of microwave radiation emitted from soils, there is a possibility 

of detecting the nearsurface water tables. 

There are some interrelated problems that must be solved in order 

to interpret measured signals correctly. These problems concern: (1) the 

physical significance of the signal emitted by a soil surface as a function 

of moisture content, (2) the surface roughness conditions that affect the 

signature, (3) the effect of vegetation medium on the soil moisture 

retrieval, and (4) the modification of the surface signal by the intervening 

atmosphere. A multi-layer incoherent radiative transfer model was developed 

by Burke et al. (1979) to analyze the 1974 aircraft experiment conducted by 

NASA in Phoenix, Arizona (Schmugge et al. 1976). The model was also 

compared with others and proved to be applicable and sufficent for a broad 

range of soil moisture conditions (Schmugge and Choudhury, 1980) The 

surface roughness effect was recently treated by Choudhury et al. (1979) 

by using a simple correction factor. It provides sufficient roughness 

correction at individual wavelengths for various observed data. One 
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shortcoming of the model is that there is a 1/ A dependence in the formu- 
lation of the correction factor; such dependence was not observed. The 
combined Burke et al. model (1979) for soil moisture (Appendix A) and 
the Choudhury et al. model (1979) for roughness correction are used 
extensively in this study for the analysis of data taken during a NASA 
conducted soil moisture experiment in 1975. Measurements of microwave 
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emissions at 21, 2.8, 1.67, 1.37 and 0,8 cm were made over furrowed bare 
fields as well as fields with various vegetation canopies. The surface 
temperature of the soil was monitored using a thermal infrared sensor. 
These instruments were carried aboard a NASA aircraft flying at an 
altitude of 500 - 800 meters over agricultural fields near Phoenix, 
Arizona. Extensive ground truth samples, in the form of moisture and 
temperature profiles were collected at the time of the aircraft over- 
passes. Section 2 summarizes the results of data analysis for both bare 
and vegetated fields. 

A comprehensive model for vegetation is presented in Section 3. Two 
perspectives are taken into account, modeling of dielectric coefficient 
for the vegetation medium and its volume scattering effect. The model 
predictions are then compared with data. In Section 4, atmospheric 
effects on soil moisture retrieval in the microwave region are summarized. 
Conclusions for this study are presented in Section 5. 


2. DATA ANALYSIS OF 1975 FLIGHT DATA 


2.1 Summary of the Experiment 

During March 1975, an aircraft mission consisting of four flights 
over the Phoenix, Arizona test site was conducted for the purpose of 
studying the use of microwave radiometers for the remote sensing of soil 
moisture. The investigators involved in this mission came from NASA, 
the Agricultural Research Service of USDA, the University of Arkansas, 
and Texas A§M University, This mission consisted of predawn and midday 
flights on 18 and 22 March 1975. There were radiometers operating in 
the wavelength range 0.8 to 21 cm. The 2.8 cm radiometer is a dual- 
polarized conically scanning radiometer operating at a fixed look angle 
of 50°. The other radiometers which were sensitive to emissions at 
wavelengths of 21, 2,8, 1,67, 1.3 and 0.8 cm were non-scanning but could 
have their nadir look angles varied. In addition to the microwave 
instruments the scientific package included a non-scanning infrared 
radiometer (10-12 urn) for measuring surface temperature. Three passes 
were taken over each field; one at a nadir angle of 0 * 0° and two at a 
nadir angle of 0 = 40° alternating the polarization sensitivity of the 
antenna. The aircraft altitude for the 0=0° pass was 800 m and for 
the 0 = 40° passes was 500 m. 

Ground measurements were made in 46 fields. Twenty-eight were 
without vegetative cover and 18 had vegetative covers of either alfalfa 
or wheat. The fields, which have an area of 16.2 hectares (400 x 400 m) , 
were arranged in pairs to provide a uniform target 800 meters wide. The 
soil moisture sampling procedures for this mission included measurements 
of the moisture content and temperature in each of the following soil 
layers: 0-1 cm, 1-2 cm, 2-5 cm, 5-9 cm, and 9-15 cm. The results 
included data from a variety of moisture conditions due to the irrigation 
and drying cycles of the fields (Schmugge, 1976). 

The data analysis was carried out for 21 cm, 2.8 cm and 1.67 cm 
measurements. .Data from 1.3 cm and 0.8 cm are not presented. An initial 
investigation showed that brightness temperatures at the shorter wave- 
lengths are relatively insensitive to soil moisture content. Responses 
to both bare and vegetated fields were studied for nadir and polarized 
conditions; 


• sc. 
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In Section 2.2, general analysis of the data are presented and compared 
with model predictions. In Section 2,3, two Stokes Parameters, P * 
h (T bv + T bh ) and Q « (T RV - T B{| ), where T gv and T BJJ are brightness 
temperatures of the vertical and horizontal polarisations respectively, 
are analyzed for optimal soil moisture sensing. In Section 2.4, the 
data are analyzed from different perspectives i'or possible retrieval of 
subs urface moisture information. 

2.2 General Analyses of Data 

2 1 cm 

Nadir and polarized responses at 21 cm for predawn (AM) and mid- 
afternoon (PM) flights are shown in Figures 2-1 to 2* 3. Average soil 
moisture contents by weight in the top 2 cm are used. For bare fields 
the correlation between soil moisture up to 20-25% and emissivity is 
high. For those fields with moisture content greater than 25%, the 
responses seem to be more variable. An explanation for this is that 
the surface roughness effect is more pronounced for moist fields. For 
example, using the estimation method developed by Choudhury et al. (1979), 
the difference in observed brighntess temperatures between medium rough 
(h * 0.3) and rough (h * 0.6) fields is only 6°K for dry fields (SM * 0%) 
but ^27°K for wet fields (SM = 25%). In general, data for bare fields 
agree well with theoretical values with an average roughness factor of 
h * 0.5. 

For vegetated fields, the response at nadir and horizontal polari- 
zation are similar. The sensitivity to soil ^aisture information is 
still strong. The general trend is that for drier fields there is a 
slight decrease in brightness response due to the vegetation as compared 
to bare fields. For moist fields, the vegetation canopy produces an 
increase in brightness temperature of about 15°K. The general increase 
in brightness temperature can be explained by the vegetation canopy 
acting as an additional thermal emitting layer. The increase is more 
for wetter fields due to their lower surface emissivities (higher surface 
reflectivity) . The increase in signature due to vegetation is further 
compensated by the smoother surfaces under the vegetation relative to 
the furrow-irrigated bare fields which would produce overall lower 
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emissivities. As a result, for dryer fields the brightness temperature 
responses over vegetated fields arc similar to, or in some cases lower 
than, those over bare fields. For vertical polarization, the vegetation 
effect is more pronounced such that the response is less sensitive to 
the background soil moisture content. The vegetation model will be dis- 
cussed in more detail in Section 4. 

2,8 cm 

The 2.8 cm data are presented in Figures 2-4 to 2-7. There are 
both vertically and horizontally polarized data at look angle of 50°. 
Three measurements were taken for each observation. The presented values 
are averages of the three measurements. In general, the correlation of 
brightness temperatures to surface moisture content in the top 2 cm for 
bare fields is better than that at 2i cm. This is due to the fact that 
at 2.8 cm, radiation comes mainly from the surface while at 21 cm, radi- 
ation from both surface and subsurface contribute to the total observed 
brightness temperature. Typically, the sensitivity of the brightness 
temperature variation to soil moisture content is about 1.5°K/% moisture 
change for data from the horizontal polarization and less than l°K/% 
moisture change for data from the vertical polarization. Data for bare 
fields also agree well with theoretical prediction with a roughness 
factor (h) ranging from 0.3 to 0.5. 

Data for vegetated fields, on the other hand, show no correlation 
to the background soil moisture. The average emissivity for the vegeta- 
tion canopy at the vertical polarization is .92 and that at the horizon- 
tal polarization is .89, both with a small variation of ±.02. Emissivity 
is defined as the ratio of the brightness temperature to the physical 
temperature using the IR observation and assuming an isothermal layer 
between the surface and vegetation. The difference between the two 
polarizations 00.03 - 0.04) can be attributed to the attenuated surface 
background which in absence of vegetation would produce a difference in 
emissivities of .07 to .1 between signatures from two polarizations. The 
observed vegetation response also agrees in general with the vegetation 
model presented in Section 4. In general, the vegetation canopy becomes 
"opaque" at this wavelength. Furthermore, brightness temperature is 
generally 20°-30°K lower than the physical temperature, a phenomenon 
that can be attributed to the scattering effect of the vegetation layer. 
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1.67 cm 


Data of 1.67 cm were analyzed in the same fashion as those of 21 cm 
and 2.8 cm. They are shown in Figures 2-8 to 2-10 for nadir, horizon- 
tally and vertically polarized observations respectively. The look 
angle for polarized observations was 40°. In general, the correlation 
of brightness temperatures to surface moisture for bare fields is 
similar to that at 2.8 cm. Radiation at this wavelength is mainly from 
the surface layer. Sensitivity of the brightness temperature to the 
surface soil moisture variation is about 1.5®K/% moisture change for 
observation from horizontal observation and less than 1°K and about 1°K 
for those from vertical and nadir observations respectively. 

Data for vegetated fields also do not show correlation to the back- 
ground soil moisture conditions. The average emissivity for the vegeta- 
tion canopy at the vertical polarization is .DO and that at the horizon- 
tal polarization is .86. The difference between two polarizations is 
similar to that at 2.8 cm. The overall lower emissivity of vegetation 
than that at 2.8 cm is resultant of the stronger scattering effect at 
1.67 cm than 2.8 cm due to the inverse wavelength dependence of the 
scattering effect which will also be discussed in Section 4. 

2.3 Polarization Effects on Soil Moisture Conditions 

The polarization information from soil moisture measurements was 
first investigated by Burke and Paris (1975) on the 1974 Phoenix data. 

Two of the Stokes parameters P * J*(Ty + T H ) and Q = (Ty - T^), where Ty 
and T u are brightness temperatures from the vertical and horizontal 
polarizations respectively, were studied to determine the polarization 
effects. It was shown that P is relatively insensitive to the surface 
roughness characteristics. Furthermore, it was also predicted that Q, 
which increases as moisture content increases, should be another useful 
parameter. The broad moisture range of the 1975 set of Phoenix data 
enabled a closer look of these effects. 

Figures 2-lla and b demonstrate the response of *$(Ty + T H ) versus 
nadir brightness temperature at 21 cm and soil moisture content respec- 
tively. As can be seen from Figure 2-lla, the nadir brightness temper- 
ature and Jj(Ty + T^) are always close to each other. This is significant 
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because it implies that in a scanning system, P may be independent of 
angle and thus provides an easy moans of analyzing data. It is also 
shown in Figure 2-1 lb that the scatter of the lj(T v ♦ Tj { ) data as a 
function of soil moisture content is loss than that of either the vertical 
or horizontal polarization. 'T’his further demonstrates the feasibility 
of utilizing P as a practical parameter for soil moisture retrieval. 

Figure 2-12 shows responses of Q * (Ty - T|j) versus soil moisture 
content at 21 cm for both bare and vegetated fields. Although there is 
a correlation between Q and soil moisture content, the scatter is too 
big for this parameter to be useful for soil moisture retrieval. How- 
ever, Q may be an indicator of surface roughness. As can be seen, Q is 
larger for vegetated fields; this is due to the smoother surfaces that 
were present under the vegetation canopy. These vegetated fields were 
not furrow-irrigated as the bare fields. Thus, it is feasible that 
utilization of the combined P and Q parameters can be quite useful for 
remotely monitoring soil conditions, 

2.4 Possible Retrieval of Subsurface Moisture Information 

The subsurface moisture information can be obtained through timely 
observations during a drying cycle from measurements from different 
wavelengths. During the 1975 Phoenix experiment, both predawn and 
afternoon data were collected on two separate days. Many of the fields 
had been irrigated about a week before the first flight so that subsur- 
face layers of these fields were moist. The second day of the experi- 
ment was conducted four days after the first, therefore allowing obser- 
vation through the drying cycle. The effective emissivities at dawn are 
plotted as a function of their values the same afternoon at wavelengths 
of 21 cm and 1.67 cm (Figure 2-13). At 1.67 cm, almost all of the fields 
have emissivities at dawn that are less than their afternoon values due 
to the evaporation process during the day. Fields with lower emissivi- 
ties both at predawn and in the afternoon can be assumed to hold water 
below the surface. Fields with the highest emissivities, which are also 
the direst, have the same values at dawn and during the afternoon, 
indicating that there is little free water available beneath the surface. 
For 21 cm, most of the radiation comes from deeper in the soil. The 
measured temperature of the soil at a depth of 15 cm generally ran 5°K 
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warmer than the surface at dawn and about 10°K cooler during the after- 
noon. The surface temperatures were generally within 2°K of the 0-1 cm 
layer value. Therefore, in addition to the moisture condition, omissi- 
vity is also affected by the relative change in the physical temperatures. 

An e Daw»* r e Aftornoon + * 05 ^ ne * 5 drawn t0 roughly allow for subsurface 
temperature effects. It is seen that for drier fields, the afternoon 
emissivity is relatively lower due to the change of surfnee-to-subsurface 
temperature gradient from predawn to the same afternoon. For fields with 
moist subsurfaces, the emissivity variation is further 'Impendent upon the 
surface moisture condition affected by the daytime evapotranspiration 
processes. 

Another way of looking at the some effect is shown in Figure 2-14 
where daw and afternoon emissivities at 1.67 cm and 21 cm are plotted 
together as a function of the surface moisture content (0»3 cm), 
Emissivities at 1.67 cm for both dawn and afternoon data all lie close 
to the same regression line. At 21 cm, the regression line is followed 
quite well for moistures greater than 15%, For fields with surface 
moisture less than 15%, the dawn observations lie above the regression 
line and the afternoon observations below. The magnitude of fluctua- 
tions of emissivity from dawn to afternoon can be used as an indicator 
of the subsurface moisture content. Thus, it is feasible that informa- 
tion of the moisture gradients of soil can be obtained from a time 
sequence of combined measurements at different wavelengths. 
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Figure 2-3 Observed brightness temperature of vertical polarization on March 18 and 22 (1975 Phoenix 
Experiment) as a function of average soxl moisture content in the top 2 cm layer for 
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of average soil moisture content in the top 2 cm layer 
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Figure 2-7 Observed vertical brightness temperature on March 18 and 22 AM (1975 Phoenix 

Experiment) as a function of average soil moisture content in the top 2 cm layer 
for X = 2.8 cm and 0 = 50° 
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Figure 2-lla The response of ^(Ty + T^) at A = 21 cm, 0 = 40° 

as a function of the corresponding nadir brightness 
temperature 
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Figure 2-llb The response of *s(Ty + T^) at A * 21 cm, 0 = 40° as a 

function of soil moisture content in the top 2 cm layer 
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Figure 2-12 The response of (T v - T H ) at X = 21 cm, 0 = 40° as a function 
of soil moisture content in the top 2 cm layer 


2-19 




EMISSIVITY (DAWN) 



0.7 0.8 0.9 1.0 


EMISSIVITY (AFTERNOON) 

Figure 2-13 Effective emissivities at A * 1.67 cm and 21 cm at dawn as a 
function of its value the same afternoon 
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Figure 2-14 Emissivities for \ * 1.67 cm and 21 cm at dawn and afternoon as 

a function of average soil moisture content in the top 2 cm layer 
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3. VEGETATION MODEL 


3.1 General Characteristics of Vegetation Canopy 

at Microwave Wavelengths 

One major purpose of soil moisture research is to predict and moni- 
tor the availability of soil moisture to agriculture application and to 
assist in the area of crop forecast. At microwave wavelengths, the 
response of .rightness temperature to surface soil moisture of bare 
fields has been well established. However, once vegetation is present, 
the signatures are contaminated to various degrees according to the type 
and coverage of the vegetation, the background condition and the wave- 
length of observation. Generally speaking, the microwave response 
observed through a vegetation canopy is as follows: 

1) the capability of radiation propagating through a 
vegetation layer is proportional to the wavelength. 

At longer wavelengths (L Band), the surface vari- 
ations can still be observed through a vegetation 
layer. At shorter wavelengths, the background sur- 
face is obscured even with small amounts of vege- 
tation cover; 

2) under conditions when surface can be observed through 
vegetation cover, the higher the background reflec- 
tivity (lower emissivity) the more "obscuration 
effect" of vegetation is expected; 

3) the amount of absorption thro! h a vegetation layer 
is directly proportional to the volume water content 
of the vegetation. The volume water content is a 
function of the plant height, density (sparse or 
cluttered) and condition (dry or wilting). For a 
forest coverage, the vegetation effect is maximum. 

Thus, the emissivity can be assumed to be close to 
one (brightness temperature * physical temperature); 

4) generally speaking, at longer wavelengths the volume 
scattering property of a vegetation medium is minimal. 
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Therefore, an approximation treating the vegetation 
medium as an additional dielectric slab above the 
surface is sufficient. This can be applied using 
the multilayer radiative transfer model developed 
for soil moisture as shown in Appendix A. The vege- 
tation medium is simply an additional layer with the 
defined dielectric coefficient. In this approach, 
dielectric properties play a major role in the 
microwave response; and 

5) at shorter wavelengths, the background surface is 
obscured in the presence of vegetation. Volume 
scattering is important and the observed brightness 
temperature is lower than the. physical temperature. 

Scattering from a vegetation layer is a function of 
the wavelength and actual structures of the plants. 

At X and K band wavelengths, a vegetation layer of 
wheat and alfalfa can produce a cooling effect of 
20-3O°K due to volume scattering such as observa- 
tions from the 1975 Phoenix experiment. Similar 
cooling effect was also observed by Kirdyashev et 
al. (1979) from fields of winter rye and corn. 

In summary, in order to model the vegetation coverage in the micro- 
wave range, two models are required along with physical properties of the 
vegetation canopy: the model of complex dielectric coefficients for 

vegetation material and the model of volume scattering for the vegetation 
medium. Modeling for the dielectric coefficient is summarized in Section 
3.2 and the volume scattering model is discussed in Section 3.3 

3.2 Models of Dielectric Coefficient for Vegetation Medium 

There are various ways of theoretically computing, or empirically 
relating, the dielectric coefficient of vegetation material at different 
wavelengths. Three representative approaches are summarized here and 
their impacts studied. 

The simplest approach is to assume that vegetation material during 
growth is mostly water (Attema et al., 1978). Using dielectric 
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coefficients for water (c wator ) and air (e air ), the equivalent complex 
dielectric coefficient of the vegetation medium ( e veg ) is simply a 
linear combination of the air-vegetation mixture: 


e veg * V v * <\ater> + V a * < E air’ ^ 

where V y and V a are volume fractions of vegetation and air respectively. 

The dielectric property of vegetation has also been estimated by 
Peake (1959) as a function of the fraction of water by weight in the 
plant (M). The simple assumption is that in the microwave range, the 
dry vegetated material has a dielectric constant of 2.5. The expression 
for the vegetated material (e v ) is thus: 

■ H * Wr * fl • *> * 2 - 5 (! - 2) 

and that for the vegetated medium is: 


veg 


V v * C* v ) 


V * 
a 


(e . ) 
v air 7 


(3-3) 


Another more comprehensive, and thus more complex, model for dielec- 
tric properties of vegetated material was undertaken by de Loor and 
Meijboom (1966). In their model, the dielectric properties of vegetated 
material was further developed as a function of both moisture by volume 

(V ) and moisture by weight (M). The expressions are: 
w 


e (Real) = 5.5 + (e - 5.5)/ (1 * f 2 x 2 ) 
v m 

e v (Im) = (e n - 5.5) * fr/(l + fV) 


(3-4) 


where ft is approximately 1.85/A (A in cm) in the microwave region and 


e ffl ■ 5 ♦ 51.5 * V w (3-5) 

V is further related to M by the densities of the solid material (ds) 

W 

and the density of water (dw) as 

M « V w /[V w ^ (1 - v w )(ds/dw)] (3-6) 

The dielectric coefficient for the vegetated medium is then derived by 
equation (3-3). 
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Table 3-1 lists dielectric coefficients at X ■ 20 cm and 3 cm 
derived from the three approaches. Assumptions are that the volume 
density of vegetation is 1%, moisture fraction by weight of the plant 
(M) is 0.8 and the density of the solid material is such that dw/ds * 4. 

TABLE 3-1 

COMPLEX DIELECTRIC COEFFICIENTS OF VEGETATION 
MEDIUM AT X = 20 AND 3 CM 



Wavelength | 

Approach 

20 cm 

3 cm 

1, Direct 

(1.8, .05) 

(1.64, ,30) 

2. Peake 

(1.64, .04) 

(1.49, .24) 

3. de Loor and 
Meijboom 

(1.30, .022) 

(1.23, .12) 


In order to show the effect of dielectric coefficient of the 
vegetation layer, the approximation of treating the vegetation medium as 
an additional dielectric slab is used. With the same assumptions for 
Table 3-1 and a vegetation layer thickness of 15 cm and background soil 
moisture content of 15%, the brightness temperature responses from three 
dielectric models are summarized in Table 3-2. 

TABLE 3t2 

BRIGHTNESS TEMPERATURES EXPECTED AT 
X s 20 CM AND 3 CM USING 
DIELECTRIC SLAB APPROACH WITH 
THREE DIFFERENT DIELECTRIC MODELS 



Wavelength I 

Dielectric Model 

20 cm 

3 cm 

1. Simple, Direct 

257. 4°K 

295. 0°K 

2. Peake 

253. 2°K 

296. 4°K 

3. de Loor and 
Meijboom 

240. 1°K 

299. 0°K 


Background Soil Moisture Content = 15% 
Vegetation Layer Thickness ■ 15 cm 
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As can bo seen, the variations in brightness temperature response 
at X * 3 cm due to different dielectric models are minimal (within 4°K). 
However, the differences at X * 20 cm are quite substantial (up to 17°K). 
Therefore, an appropriate dielectric model for the vegetation layer at 
longer wavelengths is essential. The dependence of dielectric constant 
at 20 cm is also demonstrated in Figure 3-1 which shows the similar 
variation for various vegetation thicknesses. Assumptions of plant 
density water content and biomass are the same as in Table 3-1. The 
dependence on dielectric coefficients diminishes as vegetation thickness 
exceeds ISO cm. 

3.3 Model of Volume Scattering for a Vegetation Medium 

The effects of hydrometeors on microwave radiation are closely 
related to the particle size relative to the wavelength, an effect 
frequently described in terms of the non-dimensional size parameter, 
q = 2 irr/X, where r is the radius of the particle and X is the wavelength. 
When q is very small («1), the effects are due primarily to absorption 
and can be treated as a function of liquid water content only. This 
condition applies for cloud droplets (typically of size less than 100 ym) 
at wavelengths longer than a few millimeters. For wavelengths in the 
submillimeter region, scattering effect of cloud particles has to be 
considered. For rain and surface snow, particle sizes are generally 
less than 1000 ym, therefore, volume scattering effect is important at 
wavelengths shorter than a centimeter. For a vegetation medium, it is 
expected that the volume scattering effect is important at wavelengths 
shorter than a few centimeters. The actual scattering effect of a vege- 
tation medium is yet to be determined. The reasons are (1) it is diffi- 
cult to understand the microstructure of a plant; the particles are non- 
spherical and cluttered, and (2) the distribution of plants within the 
vegetation medium is highly complex and random. There have been various 
modeling efforts of the scattering properties of various types of vege- 
tation. For example, Chuang et al. (1980) employed the concept of 
correlation lengths in the horizontal and vertical directions resulting 
in a solution dominated by the forward scattering. Fung (1979) treated 
the vegetation layer as a volume of leaves and modeled it as an inhomo- 
geneous medium with relative random permittivity function. Basharinov 
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et al. (1980) derived the vegetation model by combining surface scatter- 
ing between the vegetation and the air and the volume scattering derived 
from a stochastic model comprising spatial variations in the dielectric 
constant. 

In this study, effort was not placed in the area of defining the 
scattering properties of a vegetation layer. Instead, the resultant 
scattering effect on the observed microwave signature was investigated. 
Two aspects are considered! the response as a function of degree of 
scattering and its wavelength dependence. 

A model was developed by Burke, et al. (1979) for treating multiple 
scattering processes for microwave-infrared atmospheric retrievals. In 
this study, it is further modified for a vegetation medium. The input 
parameters include; 

1) equivalent dielectric coefficient for the vegetation 
medium and the wavelength of observation to derive 
the absorption coefficient (ya); 

2) thickness of vegetation layer to obtain the total 
absorption thickness of the medium (ya*Z); 

3) physical temperature of the vegetation medium; 

4) an equivalent single scattering albedo (w); defined 
as the ratio of scattering (ys) to extinction 

(ya + ys). The value of u> can vary between 0 (no 
scattering) and 1 (full scattering); and 

5) background soil conditions to derive the soil 
emmisivity. 

The model is summarized in Appendix B. 

Figure 3-2 is a plot of nadir brightness temperature response at 
X = 20 cm as a function of vegetation depth. The complex dielectric, 
coefficient of the vegetation layer is assumed to be (1.3, 0.022) as 
shown in Table 3-1 using the de Loor and Meijboom model. Three scatter- 
ing albedos are considered; w = 0, 0.2 and 0.5. w » 0 case is equiva- 
lent to the dielectric slab (absorption only) model described in the 
previous section. The u> = 0.2 case produces a cooling effect of the 


brightness temperature by a few degrees. The w * 0.5 case produces a 
cooling effect of , vl0 o K for vegetation thicknesses less than 50 cm and 
30° K for thicker vegetation. Comparing the data at 21 cm as presented 
in Section 2, scattering effect due to vegetation was minimal. Data 
observed from other experiments (e.g., Kirdyasher et al., 1979) also 
indicated that at L-band wavelengths the volume scattering effect can be 
neglected. 

At shorter wavelengths, the volume scattering effect is found to be 
strongly related to the scattering albedo but relatively independent of 
the vegetation layer thickness. The effect of volume scattering is shown 
in Figtre 3-3 for \ * 3 cm. The nadir brightness temperature is plotted 
as a function of the single scattering albedo (w) of the vegetation 
medium. Compared with data at 2.8 cm as presented in Section 2, a scat- 
tering effect equivalent to w = .3 to .4 was evident due to the cooling 
of 20-30° K of the observed brightness temperature from the physical 
temperature. 

Two conclusions are reached from the vegetation model and its appli- 
cation to data: (1) at L band wavelength, the background soil moisture 

information can still be retrieved; for the vegetation effect on observed 
signature, volume scattering is minimal, but the definition of the 
dielectric coefficient is important; and (2) at X and K band wavelengths, 
the background soil moisture information is generally obscured through a 
vegetation layer. Furthermore, vegetation layer also produces a strong 
multiple scattering effect within the medium, thus producing a cooling 
effect on the observed brightness temperature as compared to the actual 
physic ill t emperatur e . 
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NADIR BRIGHTNESS TEMP. (°K) 



VEGETATION DEPTH (cm) 

Figure 3-1 Response of a nadir brightness temperature at A = 20 cm 
as a function of vegetation depth for three different 
dielectric models of the vegetation layer 
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NADIR BRIGHTNESS TEMPERATURE (°K) 


280 


270 


X « 20 cm 


Ctf*0 


0.2 


260 


C<J*0.5 


250 


6$urf «0.75 

CtJ ‘SINGLE SCATTERING 
ALBEDO 


20 50 100 150 21 
VEGETATION DEPTH (cm) 

Response of a nadir brightness temperature at A = 20 cm 
as a function of vegetation depth applied to three cases 
with different scattering effect 



NADIR BRIGHTNESS TEMP. ( # K) 



SINGLE SCATTERING ALBEDO 
OF VEGETATION MEDIUM 


Figure 3-3 Response of a nadir brightness temperature at A = 3 cm 
as a function of the single scattering albedo of the 
vegetation medium 
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4. ATMOSPHERIC EFFECTS IN THE MICROWAVE REGION 


In order to define the optimum sensor wavelengths for space borne 
measurements of soil moisture , various meteorological conditions have to 
be investigated at different frequencies to test their sensitivities to 
soil moisture conditions. An atmospheric-surface Interaction radiative 
transfer model (Gaut and Reifenstein, 1971} and Burke et al . , 1979) can 
be applied to obtain the expected brightness temperature measured from 
space. The contributions of the signature include: (1) the emission of 

the attenuated background surface, (2) the upward emission from the atmo- 
sphere, and (3) the reflection of the downward atmospheric emission from 
the surface. The atmospheric parameters. contributing to the microwave 
emission signature include water vapor and oxygen. In a cloudy or rainy 
atmosphere, additional parameters of liquid water and ice crystals have 
to be considered. 

Three meteorological conditions are considered: (1) a standard 

tropical atmosphere with surface temperature of 300°K and relative 
humidity of 75%, (2) same atmosphere with a 2 km cumulus cloud from 
1.5 to 3.5 km density of 0.4 gm/m , and (3) same atmosphere with a rain 
layer up to 4 km level and rain rate of 5 mm/hr, Wavelengths of 20, 5, 

3, and 1,5 cm are tested with background surface moisture content varying 
from 5 to 35% (including a roughness factor of h * 0.4). Wavelengths 
shorter than 1.5 cm are not considered as they are not applicable for 
soil moisture monitoring due to their short skin depths. The results 
are outlined in Table 4.1. As can be seen, atmospheric effects at 20, 

5, and 3 cm are minimal under all meteorological conditions considered. 

At 1.5 cm, rain conditions can affect the sensitivity to soil moisture 
conditions. However, under cloudy conditions, the soil moisture 
information can still be obtained. In conclusion, atmospheric effects 
generally do not create ambiguities of the soil moisture information at 
wavelengths optimal for soil moisture application. 
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TABLE 4.1 


BRIGHTNESS TEMPERATURES OBSERVED FROM SPACE 
THROUGH VARIOUS METEOROLOGICAL CONDITIONS AND SOIL MOISTURE CONTENTS 


Si 

-SfiilMoi sture 
Conditions' — -- — 

5% 

15% 

25% 

35% 

a 

o 

clear 

277.0 

244,1 

213.8 

198,2 

cloud 

277.0 

244.1 

213,8 

198.2 

cs 

rain 

276.6 

243.7 

213.4 

197,8 


clear 

284.0 

244.6 

226.6 

218.8 

8 


284.1 

245.0 

227.1 

219.4 

in 

rain 

283.1 

244.5 

226.9 

219.3 


clear 

284.1 

245.1 

227.3 

219.5 

0 

cloud 

284.3 

246.2 

228.8 

221.3 


rain 

281.0 

246.3 

230.5 

223.6 

m 

clear 

285.9 

254.6 

240.3 

234.0 

IS 

cloud 

286.6 

259.3 

246.9 

241.5 

■1 

rain 

263.3 

248.1 

241.1 

238.1 
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5. CONCLUSIONS 


During tho course of this study, three tasks related to soil mois- 
ture sensing at microwave wavelengths have been undertaken: (1) analysis 

of the 1975 flight data at L, X and Ku band wavelengths from Phoenix, 
Arizona, (2) modeling of vegetation canopy at microwave wavelengths and 
(3) investigation of atmospheric effects at these wavelengths. In 
essence, information obtained from longer wavelengths (e.g., L to C band) 
can be directly applied to retrieving soil information. For shorter 
wavelengths (o.g., X, K bands) the information of surface emissivity is 
also important for atmospheric retrievals. 

Data from bare and vegetated fields included a broad range of soil 
moisture contents. Generally, tho data agree well with theoretical 
estimates from a combined multilayer radiative transfer model with 
simple roughness correction. Secondly, the polarization parameters 
P « ^(T bv ♦ T R ,p and Q * (T RV - Tg H ) are looked into. It is shown that 
P is closely related to the ridir brightness temperature and is more 
sensitive to soil moisture content than either T BV or T BH* Data of Q 
as a function of soil moisture content show . large scatter and there- 
fore may not be useful for soil moisture retrieval. Further investiga- 
tion should be carried out for studying its quantitative relationship to 
the surface roughness.. Thirdly, the possibility of detecting subsurface 
moisture information is also investigated by analyzing the predawn and 
afternoon data from combined 21 cm and 1.67 cm on two separate days 
throughout a drying cycle of the fields. 

The modeling of vegetation is carried out by considering both the 
effect of dielectric coefficient and the volume scattering characteris- 
tics of the vegetation layer. It is concluded that at longer (e.g. L, S 
and C band) wavelengths the radiation from soil can still penetrate 
through vegetation layer providing sufficient surface moisture informa- 
tion, The key factor for determination of vegetation effect is the 
definition of dielectric coefficient. Volume scattering effect on these 
long wavelengths is minimal. At shorter wavelengths, it is concluded 
that radiation from soil cannot penetrate through a vegetation canopy. 
The key element controling the signature from a vegetation layer is its 
volume scattering characteristics. 
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The Atmospheric effect at microwave wavelengths which are sensitive 
to soil moisture information is also investigated. Clear, cloudy and 
rainy conditions are considered at wavelengths of 20, 5, 3 and 1,5 cm, 

It is demonstrated that meteorological effect* car generally bo regarded 
as minimal except at shorter wavolengths under extreme weather conditions. 
There are two scientific presentations and papers accomplished under 
the contract: (1) Feasibility of Detecting Subsurface Moisture Content 

Utilizing ►fulti-Frcquency Microwave Radiometers (paper given at AGU, 
Washington, D.C., 1979), and (2) Requirements of Space-Borne Microwave 
Radiometers for Detecting Soil Moisture Contents (paper given at the 
Satellite Hydrology Symposium, Sioux Falls, S, Dakota, 1979), 

BRT is currently under another contract, as a continuation of this 
study, to evaluate statistically the surface and subsurface soil moisture 
retrievability from past data. This should provide more findings on the 
possibility of detecting subsurface moisture information and also the 
expected accuracies. 

Further research efforts should include the understanding of the 
causes of volume scattering effects of a vegetation medium by investiga- 
ting the elements that contribute to them. This could then lead to the 
possibility of identifying various vegetations and their growth condi- 
tions at microwave wavelengths. 
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APPENDIX A 

MODEL FOR MICROWAVE RESPONSE TO SOIL MOISTURE CONTENT 


Thermal microwave radiation in soils results from the random, micro- 
scopic current loops within the soil volume (Stogryn, 1970) . The inten- 
sity of radiation energy at a given point depends on the local dielectric 
coefficient and the physical temperature of the soil. Moisture produces 
a marked increase in both real and imaginary parts of the dielectric 
coefficient of soil, leading to a decrease in the soil’s emissivity. 

This effect is mainly due to a lower transmission coefficient resulting 
from an increased dielectric mismatch between the regions of radiation 
generation (soil) and the point of observation (air) . Experimental 
observations and theoretical calculations presented below indicate that 
emissivity of soils at microwave frequencies can range from >0.9 for dry 
soils to S0.& for very moist soils. 

In addition to the presence of moisture, surface roughness and vege- 
tation cover also have significant effects, generally tending to increase 
the surface emissivity. Their effects will also be discussod in Section 
3. In order to relate the microwave emissions measured by an aircraft or 
space vehicle antenna to the dielectric properties of the emitting soil, 
the following simplifying assumptions are made for the purpose of this 
paper: 

(1) the radiation within the soil is incoherent, 

(2) moisture and temperature are functions of depth only, 

(3) dielectric and thermal properties of the soil are constant 
across layers of finite thickness, and 

(4) the surface of the soil is smooth. 

A cross section of a stratified soil is shown in Figure 1. Layers 

have thicknesses AZ^ . The jth layer is bounded on the top by the jth 

surface and by the i+lth surface on the bottom. Within this layer, the 

2 2 

dispersion relation for electromagnetic wave propagation is k,. * (w/c) 

vu In this equation, w is the frequency in radians/sec, c is the 

velocity of light, p is the magnetic permeability (assumed equal to one) 

and E . = e„. + ie T . is the complex dielectric coefficient. If we write 
1 Rj v 
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k. = “(8. + ia.), the the dispersion relation gives 

K J J 

(p 2 . + p ? + 0 2 \ - (a 2 . + a ? + a 2 .) « e R . 
rxj p yj ' zj 1 i XJ yj ZJ i Rj 


20. • a. » € Ij 


and 


( 1 ) 


It can be solved using Snell’s Law to give 




( 2 ) 


At the boundary between the j and j-1 layers, radiation is partially 
reflected and transmitted. The fractions of the incident electric field 
with horizontal and vertical polarizations reflected back into the jth 
layer are given by the Fresnel coefficients. 


P H = 
P V = 


zj 


k zj~l 


k zj + k zj-l 
*j-j k zj ~ *j k zj-i 

Vl k zj + *j k zj-l 


( 3 ) 


where k . « g . + ia . . Thus 8 , a , p„ and p all depend on the complex 

ZJ ZJ ZJ Z Z M V 

dielectric coefficient and the angle 0 Q that the ray emerges from the 
soil . 

Within the first layer, the radiative transfer equation can be 
written as 



Y is the product of the density and monochromatic mass absorption coef- 
ficient. By writing the Poynting theorem in an appropriate form, it can 

be shown that y. - 2wa .(0 J/c. 

'x zi o 

I w is the intensity of radiation at frequency w, is the Planck 
emission function, In the microwave frequency range, Planck’s emission 
law reduces to the Rayleigh- Jeans equation, where is proportional to 
the temperature of the medium T. Adopting a similar scaling rule for 
I w , an effective temperature T p can be defined which is directly pro- 
portional to 1^. The subscript w is suppressed and T p refers to the 
intensity in a narrow range near o. in the p polarization state. Since 
J w is isotropic and independent of polarization, no designation is neces- 
sary. The radiative transfer equation in the first layer may be written: 


diy ± z) 



+ 


T. 


1 


(5) 


This equation can be integrated from a point just below the surface to a 
point just above the interface between the first and second layers. 
Because the dielectric properties are assumed to be constant across the 
layer 

-Y Az “ Y l Az l 

T p (l ) = Tj (1-e 1 X ) + T p (2 + )e (6) 

* 4 * 

The argument N~ implies that the measurement is made above O) or 
below (-) the Nth interface. The first term on the right hand side of 
equation 6 accounts for radiation emitted within the first layer and 
comes directly to the surface. The second term describes upwelling 
radiation at the bottom of the first layer. This in turm has two com- 
ponents: first, radiation emitted in the first layer and reflected at 

the interface between the first and second layers; and second, radiation 
transmitted from lower layers. 

Repeating the procedure of integrating the radiative transfer equa- 
tion for N layers, the brightness signature right above the surface can 
be shown as 


( 7 ) 



N / 

2?,v 

i- 1 


-*i %) 


1-c 


) 


v - -r,(« 0 ) Az i\ 

( 1+R p,i + 1 <V' ) 

■ R pj <V il j- 


An atmospheric-surface interaction radiative transfer model (Gaut 
and Reifenstein, 1971; and Burke et al, 1979) is then applied to obtain 
the expected brightness temperature measured from space, The contribu- 
tions of the signature include: (1) the emission of the attenuated 

background surface, (2) the upward emission from the atmosphere, and (3) 
the reflection of the downward atmospheric emission from the surface. 

The atmospheric parameters contributing to the microwave emission signa- 
ture include water vapor and oxygen. In a cloudy or rainy atmosphere, 
additional parameters of liquid water and ice crystals have to be con- 
sidered. 
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APPENDIX I! 


A SOLUTION TO Till: HAD I ATI VI; TRANSFIsK 
EQUATION WITH MULTIPLE SCATTERING EFFECT 


The inclusion of vegetation in the radiative transfer equation 
complicates the problem due to the addition of volume scattering effect. 
Instead of the nonlinear quadrature formula represented by equation 
(A-4), the radiative transfer equation becomes an integral equation 
incorporating the effects of multiple scattering. Numerical solutions 
to the multiple scattering equations have been available for some time 
but the procedures have been time consuming and costly to operate. With 
the variational-iterative method developed at ERT, computer times are 
short enough and costs low enough to permit the inclusion of multiple 
scattering effects in the radiance computation. To apply it to the 
microwave response to vegetation, the Planck function is replaced by the 
physical temperature and the intensity by the brightness temperature. 

Variational -Iterat ive Method 


The radiative transfer equation for a plane-parallel scattering 
atmosphere is given by (Chandrasekhar, 1960): 


P 


d !y (t , p , 4 >) 
dx 






where 




4tt 


+ 1 2 tt 

/ / P Ch . <#> ;M 1 . <#> 1 ) I (t.p 1 ,«#> 1 )dpi 1 di#> 1 

- 1 o ^ 


+ J q (t,p ,<p) 

p is the cosine of the zenith angle 
, $ is the azimuth angle 
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J 


is the source function 


is the primary oxitation source function 
p is the phase function 

and the remaining symbols were defined above. For a scattering piano- 
parallel atmosphere with the primary source emission from atmospheric 
gasses or for a plane-parallel atmosphere with an external source (the 
sun say) and isotropic scattering, the problem has axial symmetry about 
the zenith direction and the dependence on azimuth can be removed: 


dl (t,u) 
v dx 




J(t,m) 


CD 


J(T,ji) = J q (t,p) + J / X P Cm » M 1 ) I (r , p 1 ) dp 1 (2) 

where the explicit dependence on v has been suppressed. 

The phase function is normalized to represent the total energy 
scattered by a single scatter relative to the sum of the energy absorbed 
and scattered by the particle: 

1 ,2 IT 

17 LI pCMjP 1 ^ 1 )^ 1 **** ■ w < 1 (3) 

411 "1 o ° 

where w o is the single scattering albedo. For an isotropic scatterer, 
p = constant = w o . The primary exitation source function for atmospheric 
thermal emission is 

J o (x) « (1 - « 0 (t))U[T(t)] (4) 

where B[T(t)] is B^, the Planck function. In this representation of the 
radiative transfer equation, the distance from the top or bottom of the 
atmosphere is measured in units of optical depth, t, rather than distance 
or pressure. The point properties of the medium such as the single 
scattering albedo are then functions of t. 
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( 3 ) 


t(S) - / Y o (0<U 
C 


(6) 

where Y e is the volume extinction coefficient (cross section per unit 
volume) 

Y a is the absorption coefficient as defined above and C = height. 


From equations (3) , (4) and (5) , the source function can be written 


w 0 (t) 


J(t) = 1 - « (t))B(T(t)) + — j — / J(T)E 1 (|T-T*|)dx 


♦ !!2£i / 1 I(T*,M 1 )e- (T *" r)/,,1 dM 1 


where Ej is the exponential integral of the first order and 

x* is the total optical depth such that I(t*,y) is .the outgoing 
(upwelling) intensity at the lower boundary (surface) . 

The exponential integral of the nth order is defined as 


Efx) = /V^V 2 dp; n = 1, 2... 


1(t**P) has two contributions: (1) surface emission and (2) surface 

reflection. For a surface reflectivity, R, and temperature, T g , 

Ij(t*,p) = (1 - R)D(T S ) (8) 


I (t*,p) = 2*R # / I(x*,-p 1 )p 1 dp 1 


where Ifr^.-p 1 ) is the downward intensity at the surface in the direc- 
tion corresponding to p 1 . 
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Equation (13) is based on the assumption of a Lambertian surface; there 
will be an "isotropic" rofloction independent of the incident angle on 
the surface. Equation (13) can further be written as 


I (2) (T*,tO 


= 2*R*/ t J(t) 
0 


E 2 (t*- lOdt 


( 10 ) 


where E 2 is the second order exponential integral . Combining (7) , (8) 
and (9) we obtain the expression for the source function: 


(OnfT*) 

JM » Cl - » 0 (t))D(T(T)) + -Sj-i-/ JCtJBjdt-tDdt 

o 

*- £ p- ' E 2( T *- T )[(l-ll)BCT gr ) + 2.R./ T J(t)B 2 (T*-t)dt] 


( 11 ) 


The outgoing intensity at the top of the atmosphere then can be expressed 
as 


T ^ 

1(0, p) = / J(t)e" t/ji dt/p + (l-R)B(T Je" T * /M 

° T * */ 

+ 2*R •[/ J(t)E 2 (T*-t)dt3e _T * /M 
o 


( 12 ) 


The three terms on the right hand side of equation (12) are 

(1) the upward emission from the vegetation and atmosphere, 

(2) the emission of the attenuated background surface, and 

(3) the reflection of the downward atmospheric/vegetation 
emission from the surface. 

The variational-iterative (VI) approach (Sze, 1976) is used to 
solve this system of equations. The variational method depends on find- 
ing the "extremum" of a certain functional; an a priori form is used 
for the unknown function and the coefficients are found from a set of 
minimizing conditions. In essence, this method provides a direct way 
for constructing an approximate solution for the source function. The 
atmosphere is divided into subintervals and the source function is approx- 
imated as a combination of step functions in different intervals. The 
advantages of this technique are that: (1) it is fast and requires 
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little computational time to achieve satisfactory accuracy, and (2) it 
allows vertical inhomogeneity and the inclusion of surfaco reflection. 

Tho VI technique providos n diroct method for constructing an 
approximate solution to the intogral equation (11) for the sourco func- 
tion. An approximate source function can be expressed as 

J a Ct) « U a (T)Jw 0 (T) (13) 

N 

whore U (t) ■ l C.V. (t) 

a i*l 1 1 

and the V. (t) are known trial functions. The choice of trial functions 
plays an important role in the ultimate success of the variational 
method (Kourganoff, 1963). In tho variational solution employed by Sze, 
simple step functions were chosen as the trial functions. This choice 
(1) makes it simple to perform the integrals required in (11) and (2) 
the intervals can be chosen to resemble multiple cloud layers, the 
weights C i for each layer represent the average source function in 
that layer where on is the single scattering albedo for the layer. 

The total optical depth t* is divided into N-l intervals with ok 
constant over each interval. The trial function then is 

V.(T)=f 1 T j <T - T j+l (14) 

J vo otherwise 


The then are solutions of the algebraic equation 


N 

l M. .C. = f. 

fil J 1 


M. . = 6 . .At . - J 
11 1) 3 \l 


to. at. 

JL_1 


T i+1 

/ b.(t)dT 

T i 


(15) 


where 


(16) 



(17) 


°» i t 3 
*ij " ^1* i « j 


.1 *W.f 


DjCr) « / * + 1 Bj(|T-t|)dt 

m 


and f. « (^BCTCii)) + B 2 (T*-T i 3 (1-R)B(T S ) 


The variational solution is an approximation to the actual solution 
which is correct at least at one level within each layer (Szc, 1976). A 
smoothed approximation for the source function then can bo constructed: 


Jj(t) - (1-u 0 (t))B(T(t)) + e 2 (t*-t)C 1-R)B(T s ) 


,N-1 


+ ^lll l C.^TTrT D,(t) 


(18) 


j-i 

N-l 


y j 


+ w o(t) 1-L (t*-t) R l C. i|wT ( r) J ^ + 1 Ii 0 (x*-t)dt 
1 j=l J 3 T. * 

J 1 

Since the smoothed approximation is a summation over layers with oscil- 
lating residual errors, it provides a roasonablo first estimate of the 
true source function (Burke and Szc, 1977) . Improved accuracy may be 
obtained by further iterations of the integral, equation for the soui’ce 
function as: 

Vt) 


J n+1 (T) " (I-M’O)B(T(t)) ♦ -^B 2 (t*-t)( 1-R)B(T s ) 

* Mil f Bl (|t-t|)j n (t)dt 

o 

+ w o(x) D 2 (t*-t) R J T *E 2 (T*-t)J n (t)dt 


(19) 


The residue of the nth iteration is defined as 


\ ■ 


J n-1 


( 20 ) 


By specifying the maximum residue allowed, the iteration process then 
brings the source function to desired accuracy. The number of iterations 
for the system also depends on the choice of the number of step functions 
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